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Abstract. The mass radius relationship of white dwarfs, near the 
Chandrasekhar Limit, is derived for a toy model of uniform density, 
using the variational principle. A power law scaling, reminiscent of 
those found in 2nd order phase transitions, is obtained. The derived 
exponent is shown to explain the relationship obtained by numerically 
integrating the TOV equations with the equation of state for a relativis- 
tic Fermi gas of electrons. 
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1. Introduction 

A white dwarf can be sufficiently described by special relativistic physics without 
taking account of general relativity as the internal energy and pressure are both very 
much less than the rest mass density. Following Weinberg (1972) we use simple 
hydrostatics to get the equations of stellar equilibrium. 



dP GM(r)p(r) 



. P(r = 0)=P c ; (1) 
dr r z 

= +4vrrVr) , Af(r = 0) = 0, (2) 

dr 

Thus, given an equation of state p = p(p) we can obtain the mass radius relation- 
ship of the given object by integrating up to the surface indicated by p = 0. 



2. Objectives 

Assuming cold, spherically symmetric, non-rotating white dwarfs in hydrostatic 
equilibrium, the equation of state is provided by the pressure of the Fermi gas of 
electrons. The equations can be integrated analytically if the equations of state 
are polytropic, so the non-relativistic limit gives us a mass radius relationship of 
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R ~ M -1 / 3 , while the relativistic limit gives a mass, the Chandrasekhar Limit, 
independent of the radius. However the nature of the mass radius relationship in the 
vicinity of the critical mass is not revealed in the polytropic analysis. The objective 
of this work is to evaluate the radius, of a toy model proposed by Jackson et al 
(2005) near its critical mass. Its radius is shown to have a power law dependence 
on the mass. This scaling relation is used to estimate the mass radius relationship 
of a white dwarf near the Chandrasekhar Limit. The critical exponent evaluated in 
the present work is shown to fit very well with simulations for white dwarfs with a 
relativistic Fermi gas. 



3. Uniform Density Model 

Jackson et al (2005) define an uniform density star by 

p(r) = {P = 3M/AttR 3 , if r < R ; , if r > R , (3) 

and use elegant scaling arguments to transform to dimensionless coordinates de- 
fined by 

M = M/M and R = R/R . (4) 

where Mq = 2.650 M Q and Rq = 8 623 km, so that the equations of stellar struc- 
ture become simply, 

-^- = f(f;x F ,M), x F (r = 0) = x Fc ; (5) 

^r = g(r;x F ,M) , M(r = 0) = , (6) 

where the two functions on the right-hand side of the equations (/ and g) are given 
by 



3? ' 



f(r;x F ,M) = --^l -; (7) 



g(f;x F ,M) = +3r 2 4 . (8) 

Where x F = hk F c/mc 2 is the ration of the kinetic energy and rest mass of the 
electrons. Taking advantage of the scaling relations, they express the energy per 
electron in units of the electron rest energy as 

ftotal(M, x F ) = -M 2/3 x F + 3vr 2 ^ . (9) 

x F 

Where the first term is the contribution from gravitation and the second term may 
be evaluated from Fermi-Dirac statistics as 



%) = ph (1+24)^1+4 

-ln(x F + ^/T+4)] • (10) 
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Hence, the mass radius relation of the star is obtained by demanding hydrostatic 
equilibrium 

f df total (M,x F ) 
I 8xf 




4. Critical Exponent 

To evaluate a tentative value of the Chandrasekhar Limit, let us evaluate the equa- 
tion of hydrostatic equilibrium in the limit of x F » 1. 

f df total (M,x F ) \ = _M^ + 3n 2(0£{x^\ 

V 8X ? )m \ dXF X F Jm 

Substituting £(x F ), neglecting logarithmic terms with respect to linear terms and 
binomially expanding the square root up to terms of order x~2, we have 

By the choice of the dimensionless coordinates, cumbersome constants in front of 
equations have been set to 1. So, substituting 

4 = 3 <14) 

F R 3 



we get 



R 2 = M 2/3 - ^M 4/3 (15) 



This tends to zero as M tends to M cr = (3/4) 3 / 2 . Taylor expanding about this 
point, we have 



R 



1 



M 



M r 



1/2 



(16) 



It can be immediately seen from the above relation that the order parameter R has 
a critical exponent (3 = 1/2. 



5. Numerical Analysis 

The equations of structure are now integrated numerically, using a 10 6 points adap- 
tive mesh, from the central region until the pressure falls off to zero. This process 
is carried on for different central densities in the region 10 < x F < 100 to trace 
out the mass radius relationship near the Chandrasekhar Limit. 
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Figure 1. Mass radius relationship 



6. Results 

Using (3 = 1/2, the radius is written in terms of the mass as 

M 



R — R, 



constant 



1/2 

(17) 



It is best-fitted with the simulated mass radius data to give a best fit of R CO nstant = 
10149 ± 25 kms and M cr = 1.43544 ± 0.00004 Solar masses, for a Helium white 
dwarf. Investigating the parameter space, the critical exponent [3 is constrained to 
0.499 ± 0.005, which matches well with the prediction. Hence Eq. 17 gives us the 
mass radius relationship of white dwarfs near the Chandrasekhar Limit. 
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